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Institute of Physics, Jan Kochanowski University, ul. Swi^tokrzyska 15, 25-406 Kielce, Poland 

We discuss a diagonal representation of a reduced density matrix determined within the framework 
of the complex scaling method. We also discuss a possible measure of bipartite correlations in 
quantum resonance states. As an example, we consider a one-dimensional system of two bosons with 
a contact interaction subjected to an open potential well. The correlation properties of the lowest- 
energy resonance state of the system are explored over a wide range of the inter-boson interaction 
strength, including the Tonks-Girardeau regime. 


I. INTRODUCTION 


eluding remarks. 


Quantum entanglement is a fundamental feature of the 
quantum world and has attracted significant research at¬ 
tention. New experimental techniques have opened up 
opportunities for practical applications of quantum en¬ 
tanglement in various branches of information technol¬ 
ogy [Ei- The fields of quantum teleportation, quan¬ 
tum cryptography and quantum computation have grown 
particularly rapidly in the last few years. Entanglement 
is also used as an alternative measure of correlation in 
systems of interacting particles Considerable efforts 
have been made to understand quantum correlations in 
bound states of model systems such as the Moshinsky 
atom quantum dot systems [Ml or ultra-cold 

boson systems (T3 - [ll . Moreover, in recent years, the 
helium atom and helium-like ions have been extensively 
studied in this context [iMl . For an overview of the re¬ 
cent developments in studies of entanglement in quantum 
composite systems, see [H. 

However, relatively few attempts have been made to 
improve understanding of the correlation properties of 
systems that exhibit metastable states. Various theoret¬ 
ical methods to determine resonance energies and life¬ 
times can be found in the literature. Among these, the 
most popular are the complex scaling method (CSM) 
^1 and the real stabilization method [11. Treatment 
of the entanglement of resonance states with the CSM 
was proposed in [1 , where both a complex-scaled density 
operator and a complex linear entropy were introduced. 
Within the framework of this formalism, the resonance 
states of two-electron Gaussian quantum dots have re¬ 
cently been analysed from the perspective of quantum 
information [H. Nonetheless, as far as we know, no 
study has discussed in detail the diagonal representation 
of the complex-scaled reduced density matrix. This gap 
in the literature provides the motivation for the present 
Letter. 

The remainder of this Letter is structured as follows. 
Section im briefly outlines the CSM. Section Hill discusses 
the diagonal form of the reduced density matrix deter¬ 
mined in the framework of the CSM. Section HVl discusses 
possible correlation measures for resonance states. Sec¬ 
tion El focuses on the correlation properties of systems 
that contain two interacting bosons trapped inside an 
open potential well. Finally, Section El presents con- 


II. THE COMPLEX SCALING METHOD 

The CSM is a powerful tool for searching for the res¬ 
onant par ameters of a system that supports metastable 
states [^. The utility of this method is that the res¬ 
onant states can be treated by applying the methods 
used to compute bound states, for example, a flnite-basis- 
set approximation [s^, Issjl . Here we only briefly outline 
the complex scaling formalism. Details of the computa¬ 
tional technique of the CSM can be found in an excellent 
overview [3ll |. 

In the CSM, the Hamiltonian H = T(x)-|-V(x), where 
T and V are kinetic and potential energy operators, re¬ 
spectively, is transformed by coordinate transformation 
X !->■ xe^® into 


fye ^ g-2/e-f F(xe^®), (1) 

where 6 is the so-called scaling angle. A key aspect of 
the CSM is that the complex-scaled Hamiltonian is a 
non-Hermitian operator and the inner product is defined 
as 


■C tplip >= / ilj{x.)ip{x.)dx, 

J all space 


( 2 ) 


<C x|(/5 (/5(x) and ^ '0|x >>= V^(x). The right {R) 

and left (L) eigenstates of are defined by 
Wilijjf ^ and ^ respectively. These 

equations can be turned into algebraic problems by ex¬ 
panding ^ and <C 4’i\ as linear combinations of 
\tpk ^ and Pk\, respectively, where {pi} is a com¬ 
plete set of orthonormal functions, <C Sij. 

Thus, the first equation becomes the eigenvalue prob¬ 
lem (H - IW)R = 0 with H = [< p*\H^\pk »], the 
eigenvalues of which Wdiag = diag{Wo, Wi ,...,) (if H is 
diagonalizable) are given by Wdiag = R^^HR, where 
R is the eigenvector matrix of H. The second equation 
becomes the eigenvalue problem (H^ — IW )L = 0. Be¬ 
cause Wd.ag = fVL, = (R-'HR)^ = (HR)^(R-i)^ = 
R^H^(R"^ )^, it immediately follows that the eigen¬ 
vector matrix of the matrix is L = (R^^)^ and 



2 


Wi = W^. As a result, we obtain 

< xlV'f »='i/'*^(x) = ^(R)fci</jfc(x), (3) 

k 

and 

< »= V'j^(x) = (4) 

k 

Because R“^R = I, the family {V’f'iV’f'l forms a com¬ 
plete set of orthonormal functions with respect to the in¬ 
ner product (12), that is, -C ^ij- In particular, 

if {(^i} is a real basis and H is a symmetric matrix (H = 
H^), that is, Wdraa = = (R^HR)^ = R^H^R 

(R~^ = R”^), then the right and left wavefunctions 
can alternatively be expressed in the same form, namely 
V’f = Efe (R)fci</5i(x), where R'^R = I. However, in 
the case of = H* (H is Hermitian), that is, Wdiag = 
W*^^g = (RIRR)* = R^H*R* = R^H^R* (R-^ = 
RI), they can be written as ^/^(x) = X]fe(R)fe*‘/^fc(^) nnd 
'0f'(x) = (R*)fe*V^fc(^)i where R^R = I. Note that 

our conclusions coincide with those of [3l|. 

The resonances appear as 0-independent complex 
eigenvalues H4 with Im[H4] < 0 [^, and the reso¬ 
nance energies and lifetimes T^ are obtained as 

= Re[lTfe] and T^ = —2Im[lTfc], respectively: 

Wk = Er-I^. ( 5 ) 

III. COMPLEX-SCALED REDUCED DENSITY 
MATRICES 

Suppose we divide a system into two parts, A and B. 
Let us express the right and left wavefunctions of a given 
resonance state as follows: 

1^^ >= X] (6) 

b 

and 

b*\- (f) 

kl 

Here, <C » = 1, where {uj} and {bj} are bases 

of square integrable orthonormal functions for the sub¬ 
systems A and H, respectively, = 6ij, 

B<C = Sij. Following Q, we define the den¬ 

sity operator as 

= 1^/,^ X(8) 

The reduced density matrix of subsystem A is thus ob¬ 
tained by tracing out the B degrees of freedom, which 
gives 

PA=tTB[p'^^] = J2P^k\a^:^AA<^ all, ( 9 ) 

ik 


where e^e^y Its right and left eigenstates are 

defined by /5 a»a= »a and ^<C u^\pA = 

^<C u^|A„, respectively. In strict analogy with Section 
El we conclude that 

|u„ f^A= ^ A -, (10) 

3 

and 

A«u^| =^(V-1)J„ ^«a;|, (11) 

3 

where V is the eigenvector matrix of the matrix pA = 
[a< a *\ pA\ak >a] = [pfkl with the eigenvalues (A„ = 
A„ = A;^), ^<C »A= 5nm- Following on from the 

above, pA can be expressed in diagonal form as follows: 

/5a = >aa< M^l, (12) 

n 

[a< u^IpaIm^ >a] = diag{X^,Xt, ■■■,)■ Because 
= 1, the normalization condition gives 
X)nAn = 1 or equivalently tr[pA] = 1- Note that the 
eigenvalues of pA are generally complex numbers, except 
for the case pA = Pai that is, when pA is a Hermitian ma¬ 
trix. Once again in strict analogy with SectionEl we can 
conclude that if {oi} is a real basis and pA = Pa’ then 
the right and left eigenstates of pA can be expressed in 
an identical form, namely 

A< >a=a< x\u^ »a= ^ {^)knak[x), (13) 

k 

V^V = I. Analogously, we obtain the reduced density 
matrix for subsystem B: 

/5b =trA[/5^^] =^p,^|0i»BB< 0fc|, (14) 

ik 

pfk = its diagonal form, 

PB = '^Xl^\v^ :>bb<^v^\, (15) 

n 

with 

k^»B=^(U),n|6,»B, (16) 

3 

and 

B«^nl =E(U“')Jn B«^*l, (17) 

3 

where A^ are the eigenvalues of the matrix pB = 
[b< 5*|pB|5fe >b] = [pfkl 3'nd U is the corresponding 
eigenvector matrix, B< »B=<^ nm ■ 

Let us rewrite the matrix pA = [pf^h pfk = 'Ej 
and the matrix ps = [pH], pH = Ej as pA = 


3 


and ps = (e^)^e^, respectively, where = 
It is known that if C and D are square complex 
matrices of the same size, then the matrices CD and DC 
have the same eigenvalues [s^- We thus conclude that 
Pa and p^ = (e^)^e^ have a common set of eigenvalues. 
Hence, bearing in mind that ps has the same eigenvalues 
as , we arrive at the conclusion that the eigenvalues of 
PA and pb are identical, Xf = Xf = Xi. 

IV. CORRELATION MEASURES 

As discussed by Moiseyev [3l|, the real and imaginary 
parts of the mean value of a given complex-scaled oper¬ 
ator, <C Q i^LlQl'ipR give the average value 

of the quantity under consideration, and its uncertainty, 
respectively. It is easy to see that the average value of 
any operator acting on one of the subsystems, let it be 
A, is given by 

< Qa >=< iPlIQaH^r >= Y. efeeg al\QA\ai > 

ijk 

Noting that = pfk^ obtain 

< Qa >= Y al\QA\ai >a= ti[pAQA]- (19) 

ik 


of 5, (Sun) as the entanglement entropy of a resonance 
state is problematic. Despite of this fact, we call S and 
Siin complex entropies and propose them as measures of 
correlation between the subsystems A and B. Thus, we 
identify the real and imaginary parts of S, {Sun) with the 
amount of correlation and the uncertainty of this amount, 
respectively. Note that the complex linear entropy of a 
resonance state was first introduced in Q. 

In a way that is analogous to the standard quantum 
theory, the resonance state that is factorized as a prod¬ 
uct of states can be regarded as uncorrelated, this cor¬ 
responds to the case in which only one eigenvalue is 
nonzero, Xi = 1, giving S = Sun = 0. Deviations from 
such a state can be characterized by the complex entropy^ 
as discussed above. 


V. EXAMPLE: TWO-BOSON SYSTEM 

Thus far, we have kept our discussion quite general, 
an example, we now consider a simple model system 
composed of two identical bosons interacting via a con¬ 
tact potential of strength g. For the sake of simplicity, 
we model an external potential by v(x) = 0.5x‘^e~^ so 
that the system does not exhibit any bound state. The 
resonance parameters of this system can be found from 
the following complex-scaled Hamiltonian 


Let us now address the question of how to characterize 
the correlation in resonance states. Generally, according 
to the standard quantum theory, the von Neumann (vN) 
entropy S = — tr[p^ b^apa,b\ (38| and the linear entropy 
Siin = 1 — tr[/9^ are used to quantify the degree of 

entanglement in composite quantum systems. In a strict 
mathematical sense, the values of S and Sun can be given 
by S' = -(InpA.s) and Sh« = (1 - Pa,b) = 1 - {pa,b) 
(^ . where (...) is the conventional inner product and i 
is the identity operator from the appropriate basis set. 
Accordingly, in the resonance case, we have 

S = - < InpA.B »= -ir[pA,B^npA,B\ = 

= -^AdnA, (20) 

i 

and 

Siin = 1 - < PA,B >= 1 - ^^[P%b\ = 

= l-^A^ (21) 

i 

Thus, the real and imaginary parts of S, (Sun) 
can be identified as the mean value of the operator: 
—\n.pA,B, (1 — PA,b) and its uncertainty, respectively. In 
particular, the real part of ^ pa,b gives the 

average of the probability Ai, whereas the imaginary part 
describes its uncertainty. However, because the entropy 
cannot rigorously be treated as the average value of a 
quantal observable [4l| , the interpretation of the real part 


2 —219 

i=l * 

Here we apply the basis function method and diagonalize 
the matrix Hamiltonian H® = ((j)nm\H^\(j)ij) in a basis of 
permanents constructed from the one-particle orthonor¬ 
mal basis set, 


= Sij[ipi{xi)ipj{x2) + 'tlJj{xi)'tlJi{x2)],{i > j), (23) 

where sa = 1/2 and = 2“^/^ for i ^ j. We choose 
as the one-particle basis the wave functions of a simple 
harmonic oscillator. 


tpi{x) 



(24) 


so that the basis (l23l) is real. Hence, and because H® is 
symmetric, then, in accordance with Section [Hi the right 
and left wavefunctions of a given resonance state can be 
written in the same form, 


X^’'”(a;i,a;2) = ^ (ii, 0:2) = E(^)b'^i(^i)V'j(^ 2 ), 

i>j ij 

(25) 

with where {ry} is the corresponding 

eigenvector of H®, and (e)^^ = ru and (e)y = 2“^/^ry 
(2“^/^rji) for i > j {i < j). From here on we denote 
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jg gg^gy ggg^ ^]^g reduced density ma¬ 
trix for particles 1 or 2 is 

P12 = X '4}j\. (26) 

Because the matrix e is symmetric, its eigenvector matrix 
V (V^ = ) and eigenvalues D = diag{do,di,...,) 

satisfy e = VDV^. Hence it is easy to infer that = 
VD^V^. After rewriting these formulas as 

(e)., = (27) 

n 

and 

=J2{VhndliV),^, (28) 

n 

substituting them into Eqs. (ESI) and Ei), respectively, 
and performing some straightforward algebra, we arrive 
at 

X (^11^2) — ^ ^ ^n^n(^l)^n(^ 2 )? (^^) 

n 

and 

pi 2 (x,x ) =< a;|/3i2|x A„u„(x)u„(x ), (30) 

n 

where 

Unix) = y] (V)knipk{x), (31) 

k 

and Xji ^nl^m ^nm' 

The resonance parameters are determined at the angle 
d = Oopt at which the eigenvalues of exhibit the most 
stabilized characters with respect to 6. We recall that 
the resonance energy and lifetime F are obtained 

from the stable eigenvalue W as = Re[bF] and F = 
—2Im[lF], respectively. We find that the set of basis 
functions (ESI) constructed from the 90 lowest one-particle 
orbitals (I^Sl) is sufficiently large to obtain a good estimate 
of the parameters of the lowest-energy resonance state, 
at least over the range of (/ = 0 to g = 45. Moreover, 
in this range the optimal value of the parameter 9 is 
approximately Oopt = 0.2, regardless of g. 

We now examine the properties of the system. Let 
us first briefly discuss the special cases g = 0 and 
g —>• 00 , which correspond to the non-interacting case 
and the Tonks-Girardeau (TG) regime [i^, respectively. 
In these limiting situations, the positions of the low¬ 
est resonance states are given by = 2Wo and 

_j_ respectively, where Hb, Wi are the 
lowest-energy resonance positions of the corresponding 
one-particle system, which we find numerically to be at 
Wo a: 0.411 - 0.00267, Wi « F014 - 0.1257. Here we use 
the complex linear entropy Sun as a measure of the corre¬ 
lation, Siin = 1 — trp^ 2 - Figs. [T] and E] show our numer- 



g 

FIG. 1: The resonance energy and lifetime as functions of g. 



g 


FIG. 2: The real and imaginary part of Sun as functions of g. 


ical results for W and for Sun as functions of g, respec¬ 
tively. The horizonal lines in both figures indicate the 
results for the TG system, Ri 1.425, Ftg ~ 0.254, 
^Hn ~ 0.34 — 0.047, where the last result was determined 
from the resonance TG wavefunction constructed as 

XTG(a^i,a^2) = sgn[x2 - (xj), (32) 

where ^g®^ and resonance or¬ 

bitals of the one-particle system corresponding to Wq ~ 
0.411 — 0.00267 and Wi ~ 1.014 — 0.1257, respectively, 
<C In the non-interacting case, we 

have £;^®b ~ 0.822, Fs=o r; 0.0104 and Xg=1](^T^2) = 

4 ‘o‘^^ixi) 4 >o^^ix 2 ), which gives Sf^~° = 0, reflecting the 
fact that there is no correlation between the particles. 
We can observe how the results obtained for finite val¬ 
ues of g converge to those for the TG system as g is 
increased, which, in particular, confirms the correctness 
of our calculations. In fact, the system starts to exhibit 
the behaviour of the TG system after exceeding a value 
of g R 40. We conclude from our results that the larger 
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the value of g, the higher the correlations produced by 
two bosons, which is attributed to the fact that the real 
part of Siin increases with the increase in g. As expected, 
the effect of changing g becomes less pronounced as g be¬ 
comes larger and disappears in the limit g —>■ oo. 

VI. CONCLUDING REMARKS 

We have discussed in detail the diagonal representa¬ 
tion of a reduced density matrix determined under the 
framework of the CSM. Moreover, we discussed the quan¬ 


tification of bipartite correlations in quantum resonance 
states by means of the complex entropy. We also con¬ 
ducted a comprehensive study of the lowest-energy reso¬ 
nance state of two interacting bosons trapped inside an 
open potential well. Among other findings, our results 
show the dependence of the complex linear entropy on 
the inter-boson interaction strength g. Its real and imag¬ 
inary parts have monotonically increasing behaviours as 
g increases and tend to constant values in the TG limit. 

We hope our study will stimulate broader discussions 
of correlation in quantum resonance states. 
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